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ABSTRACT 

The first hght from a supernova (SN) emerges once the SN shock breaks out 
of the stellar surface. The first light, typically a UV or X-ray flash, is followed 
by a broken power-law decay of the luminosity generated by radiation that leaks 
out of the expanding gas sphere. Motivated by recent detection of emission from 
very early stages of several SNe, we revisit the theory of shock breakout and 
the following emission, paying special attention to the photon-gas coupling and 
deviations from thermal equilibrium. We derive simple analytic light curves of 
SNe from various progenitors at early times. We find that for more compact 
progenitors, white dwarfs, Wolf-Rayet stars (WRs) and possibly more energetic 
blue-supergiant explosions, the observed radiation is out of thermal equilibrium 
at the breakout, during the planar phase (i.e., before the expanding gas dou- 
bles its radius), and during the early spherical phase. Therefore, during these 
phases we predict significantly higher temperatures than previous analysis that 
assumed equilibrium. When thermal equilibrium prevails, we find the location of 
the thermalization depth and its temporal evolution. Our results are useful for 
interpretation of early SN light curves. Some examples are: (i) Red supergiant 
SNe have an early bright peak in optical and UV flux, less than an hour after 
breakout. It is followed by a minimum at the end of the planar phase (about 10 
hr), before it peaks again once the temperature drops to the observed frequency 
range. In contrast WRs show only the latter peak in optical and UV. (ii) Bright 
X-ray flares are expected from all core-collapse SNe types, (iii) The light curve 
and spectrum of the initial breakout pulse holds information on the explosion 
geometry and progenitor wind opacity. Its spectrum in more compact progen- 
itors shows a (non-thermal) power-law and its light curve may reveal both the 
breakout diffusion time and the progenitor radius. 
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Introduction 



A breakout of a shock through the stellar surfac e is predicted to be the first electro- 



magnetic signal heralding the birth of a supernova rSN) JColgatelll974j:lFalklll978l 



1992 



Klein fc Chevalier 



Matzner fc McKee 



1978 : llmshennik. Nadezhin fc Utrobiru Il981t lEnsman fc Burrows 
19991 ). Before breakout the shock is propagating through the opaque stellar envelope. The 
shock is radiation dominated (i.e., the energy density behind the shock is dominated by 
radiation) and it accelerates while propagating through the decreasing density profile of the 
envelope, leaving behind the shock an expanding radiation dominated gas. Following the 
shock breakout, photons continue to diffuse out of the ex panding stellar envelope producing a 



long l asting emission that slowly decays with time (e.g . . iGrassberg. Imshennik fc Nadyozhin 
197ll : IChevaherlEgTei Il992l : Ichevalier fc Franssonll2008l : IPiro. Chang fc Weinberg|l2009[ ) . The 
typical frequency of the breakout emission ranges from far ultra-violet to soft 7-rays in core 
collapse SNe, and as we show here, is in 7-rays in type la SNe. The typical frequency of the 
following emission decreases to the visible-near UV bands after a day. The energy released 
during the breakout increases with the progenitor radius and can reach ~ 0.1% of the SN 
explosion energy in a red supergiant. The luminosity of core collapse SNe after a day is 
~ 10^^ — 10^^ erg/s. Thus, the shock breakout and the emission through the first day can 
be detected out to the nearby Universe, but without any preceding knowledge of where to 
look, their detection is challenging. Nevertheless, the search worth the effort as this emission 
bears direct information on the properties of the progenitor and the explosion, which are 
difficult to obtain in any other way. 

The development during the recent decade of sensitive UV, X-ray and soft gamma-ray 
detectors, with relatively large fields of view, lead to the discovery of several shock breakout 
cand i dates (ICampaiia et al. II2006I : ISoderberg et al. iGezari et al. ll2008l:ISchawinski et al. 



2008 : Modjaz et al. 20091 ). Motivated by these, and by the rising potential for future detec- 
tion of shock breakouts from various progenitors, we revisit this topic. We develop an analytic 
model that provides light curves (luminosity and temperature) starting from the breakout, 
through the quasi-planar expansion phase to the spherical expansion phase, until recombina- 
tion and/or radioactive decay start playing a significant role. These phases were explored in 
previo us works, where t he most updated analy t ic stu dy of the spherical phase was carried- 
out by IChevalierl (Il992[): IChevalier fc FranssonI ( 120081 ) , and IWaxman. Meszaros fc Campana 
( I2OO7I ): iRabinak &: WaxmanI (120101). The study of the planar phase was carried out only 
very recently by iPiro. Chang &: Weinbergl ( 120091 ) in the context of Type la shock breakout. 
The advantage of our model is that we follow the photon-gas coupling within the expanding 
gas. At each stage of the evolution we find the location at which the observed temperature 
is determined. We find whether the radiation at this place is in thermal equilibrium or not 
and calculate the observed temperature. It turns out that the temperature evolution during 
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the planar phase and the early spherical phase depends strongly on the thermal equilibrium 
of the radiation just behind the shock in the breakout layer. In radiation don iinated shocks , 
radiation is out of thermal equilibrium if the shock velocity is high enough (jWeaverlll976l ). 
which is th e case in shock breakout from more compact progenitors and more energetic 
explosions (IKatz. Budnik &: Waxmaru 120101 ) . We provide the first calculation (analytic or 
numerical) of the light curve in case that the observed radiation is out of thermal equilib- 
rium at the source. We also carry out the first analytic calculation of the evolution of the 
location of the thermalization depth, through the different phases, when the radiation is in 
thermal equilibrium. We use our model to examine the light curve and spectrum of the 
initial pulse that is strongly affected by light travel time effects and opacity of the progenitor 
stellar wind. Finally, we use our model to explore the properties of early SNe light curves 
resulting from various progenitor types including red supergiants (RSG), blue supergiants 
(BSG), Wolf-Rayet stars (WR) and white dwarfs (WD). We provide simple formula of early 
SN light curves for these different progenitors. 

We present our model and its general results in section |2l The bolometric luminosity 
and spectrum of the initial pulse are discussed in section El Early SNe light curves resulting 
from various progenitors are presented in section m A reader that is interested only at the 
final light curves should go directly to this section. In section |5] we compare our calculations 
to previous analytic and numerical studies. We summarize our main results in section [61 We 
ignore in this paper any cosmological redshift effects. 



2. Luminosity and observed spectrum following shock breakout 

SN explosion drives a radiation dominated shock (i.e., the internal energy in the shock 
downstream is dominated by radiation) that propagates and accelerates through the decreas- 
ing density profile of the stellar envelope. The pre-explosion density profile near the stelar 
radius, i?*, can be approximated by a power- law p oc (i?* — r)" where r is the distance from 
the star center and n ^ 1.5 — 3, depending on the stellar properties. In such medium the 
blast-wave assumes a self similar profile in which the shock velocity v oc p^^, where /i is 
a very weak functi on of n and f or the relevant range of n can be well approximated as a 



constant /i ~ 0.19 (jSakurailll960l ). This purely hydrodynamic solution, which neglects heat 
conduction, holds as long as the distance to the stellar edge, —r, is larger than the width 
of the radiation dominated shock, or equivalently as long as the optical depth for photons 
to escape the stellar envelope, T(r), is larger than c/v{r), where c is the light speed. Once 
T c/v the shock "breaks-out", leaving behind an expanding hot gas-radiation sphere. We 
denote the conditions at the front shell of gas at the breakout time (where r ~ c/f at the 
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breakout) by subscript o, and call this shell the breakout shell. Given above assumptions 
our problem is completely defined by the following parameters of the breakout shell at the 
breakout time. The parameters of the inner gas layers are directly related to it. 

• Vq: the velocity of the shock when it breaks out. 

• mo: the mass of the breakout shell. 

• do: the initial width of the breakout shell. 

• R^: the radius of the star. 

• n: the power law index describing the pre-explosion density profile, typically 1.5 < 
n < 3 

Other characteristics of the breakout shell can be calculated from the above. For example, 

• tq = Kj-mo/i?^ is the initial optical depth for Thompson scattering, where is the 
Thompson cross section per unit of mass. 

• to = do/vQ is the initial dynamical time of the breakout shell. It is also its initial 
diffusion time. 

• Eq = mof is the initial thermal energy of the breakout shell. 

Conceptually, it is useful to treat the whole expanding envelope as a series of successive 
shells. Any sh ell deeper than the breakout shell, can be characterized by its mass m > mo. 



Sakurail (Il960l )'s hydrodynamic solution discussed above indicates that for these deeper shells 



• V = fo(m/mo) "+i is the velocity of the shell of mass m. 

• di = dQ{m/mQ)~ is the width of a shell of mass m at shock breakout. 

• Ei = mv = Eolm/mo) "+i is the internal energy of a shell of mass m at shock 
breakout. 

The subscript i, for time dependent quantities, indicates initial values. 

Initially, after shock crossing but before significant expansion, the thermal and kinetic 
energies in each shell are equal, hence the expression for the initial internal energy given 
above. Following shock breakout the stellar envelope expands and the width of each shell 
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increases as it accelerates on the expense of the gas thermal energy. The acceleration ends 
once the shell width is mu ltiplied several tinies, ap proaching a coasting velocity which is 
about twice its initial one (IMatzner fc McKeelll999l ). Hence, up to a factor of order unity, 



the initial velocity right after the shock equals its final coasting velocity, v. The internal 
energy on the other hand falls as the volume of the shell increases. 

Following shock breakout, photons continue to diffuse out of the expanding stellar en- 
velope. Below we derive the luminosity and typical frequency of the diffusing photons. We 
focus on the early phases, where the gas ionization fraction is high and the opacity is dom- 
inated by Thomson scattering (e.g. over free-free absorption). This assumption breaks as 
soon as the observed temperature falls below about 1 eV (the gas density at this point is 
XO~io — 10~^ g/cm^). We also neglect the drop in the number of free electrons due to re- 
combination of fully ionized He (which takes place at a few eV), and assume that the energy 
injection due to recombination and radioactive decay are negligible. These conditions prevail 
from the shock breakout up until about a day after the explosion. 

First, we derive the observed bolometric luminosity as a function of time. It is dic- 
tated by the hydrodynamics and Thompson scattering and does not depend on the thermal 
coupling between the diffusing radiation and the gas (at the stages of interest the shells 
are coasting and the hydrodynamics is decoupled from the radiation). Then, we discuss the 
radiation-gas coupling and derive the observed spectrum, both when the diffusing radiation is 
in thermal equilibrium and when it is not. We do not provide the radius of the photosphere 
(i.e., r = 1) since it has no observable signature in the specific system at hand. Instead 
we provide the radius and mass of the location at which the luminosity is determined (see 
definition below), r and m and the radius and mass of the location at which the observed 
temperature (color) is determined, and rrici. 



2.1. Luminosity 

At any given time the observed luminosity is dominated by the innermost shell out of 
which photons can effectively diffuse, i.e. where the diffusion time is comparable to the 
time past since breakout, which is also the dynamical time. The reason is that the diffusion 
time from inner shells is longer and therefore their photons are confined and cannot be seen, 
while the radiation from outer shells already escaped the expanding gas and were observed 
at earlier times. The criterion that the diffusion time is comparable to the dynamical time 
is r ^ c/v. Thus, the bolometric luminosity at time t is determined by the radiation energy 
in the shell that satisfies r(t) ~ c/v. We call this shell the luminosity shell, and denote x as 
the parameter x of this shell. The luminosity is given then by E /t. 
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To follow the evolution, we measure the time t from breakout. The relevant parameters 
of a shell are the width d, velocity v, internal energy E, optical depth r, and photon diffusion 
time td, all of which, except for v, are time dependent (f only increases by a modest factor 
from its initial value, and is approximated here as constant). The following relations holds 
for each shells] at all times t > to: 

m 

T = Kt — 

Td 
= — 



d = di + vt ^ vt 

r ^ R^Wt (1) 

Our equation for the internal energy above assumes adiabatic expansion of a radiation dom- 
inated gas. It is therefore not valid for outer shells where t > td which cooled radiatively. It 
is applicable from the luminosity shell and inward. 

The hydrodynamic evolution has two phases - a planar phase and a spherical phase. 
The evolution of a shell is approximately planar as long as its radius did not double and it 
is spherical at later times. As we show below the breakout shell is also the luminosity shell 
while the evolution of the breakout shell is planar. Thus we separate the temporal evolution 
into two asymptotic regimes, (to <)t ^ tg and t ^ tg, where ts = R*/vq is the time of 
transition between the planar and spherical geometries of the breakout shell. 

During the planar phase r ^ R^, implying a constant r. Since d oc t we obtain that td/t 
is constant in time for all the shells during their planar evolution. As we have seen above 
at the beginning of the planar phase (i.e., t = tg), the diffusion time equals the dynamical 
time only for the breakout shell. Therefore, the breakout shell is also the luminosity shell 
(i.e., m ^ mo) throughout the planar phase. The implication is that throughout the planar 
phase we continue to observe photons from roughly the breakout shell. Adiabatic cooling 
then dictates 

E{t < ts) ^ Eo{t/to)-'/\ (2) 

1.38n + l 

During the spherical phase r vt and the opacity is, t ~ (^) ^ = tq "^^ (yj- 

The luminosity shell satisfies, by definition, td = t, which is equivalent to r = c/v. Com- 
bining these two equations for r, we find that at t > ts the luminosity shell evolves with 



^The last approximation d ^ vt holds at only for mp. For other shells it holds dX t > di/v. However, 
these shells are irrelevant before that time so this approximation is valid for our analysis 
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^ 2(n+l) ^ 0.38w ^ O.Slw+l '> 2 ^ 0.38n 

m oc i i i9"+i , V (xt i i9"+i , r oc f i wn+i ^ c^. oc i i isn+i and r oc ^i isi+i . Using these relations, 
and following the adiabatic cooling of the shells we obtain 



-1/3 



1.29ra+5 



+ \ 3(1.19n+l) 

E{t>Q^Eo(f] i^J . (3) 



The observed bolometric luminosity, given by E/t, is therefore: 

-4/3 




Lobs ~ \ . s -4/3 . s 2-28n-2 (4) 

i 1 / ^ "'/•^ / ^ \ 3(1.19n-H) 



t < ts 
t > ts 



The ratio Eo/to is simply the initial luminosity Lq. This bolometric luminosity falls as t~^^^ 
in the planar phase and falls much more modestly as to t~^-^^ (for 1.5 < n < 3) during 
the spherical phase. The total energy released in a logarithmic time unit, Lt, increases once 
the spherical phase sets on. The reason is that most of the energy in the shock is left behind 
it, in deeper more massive shells, as the shock accelerates (indeed we have seen that Ei is 
an increasing function of m). Therefore inner shells, which were previously opaque, become 
transparent during the spherical phase and those contain increasing amount of energy. 



The radius and mass of the luminosity shell are: 



t < ts 



p ; o.siw-n (r.\ 

* 1 / t \ 119n-Hl , , ) \'-'J 



t < ts 
t > ts 




2.2. Spectrum 

Throughout the early phases that we consider here, the radiation dominates the heat 
capacity. Therefore, the luminosity derived above is independent of the thermal coupling 
between the photons and the gas (even a full coupling which equates the gas and radiation 
temperatures everywhere has a neghgible effect on the total radiation energy). However, 
the optical depth of the luminosity shell is always much greater than unity (r^ — tq ^ 1 
and T only increases, or remains constant, with time). Thus, the escaping photons would 
have many interactions with electrons as they cross the envelope until finally escaping from 
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the layer where r = 1. Therefore, while the total radiation energy is independent of the 
photons-electrons interaction, the typical energy of each of the escaping photons is dictated 
by their coupling to matter along their way from the luminosity shell outward. 

At any given time the photon luminosity at r > r is constant (given by equation H]). 
Therefore, there are two extreme possibilities. (1) No radiation process can significantly 
change the number of photons. Then, since photons dominate the heat capacity, both 
luminosity and photon number flux are fixed (independent of r for r > r). Thus, the 
typical photon energy must be fixed. The typical temperature is thus given by the photons 
temperature T of the luminosity shell. That temperature, in turn, would be dictated by 
the limited ability of radiation processes in the luminosity shell to create photons. (2) Some 
radiation process can create new photons, which can share the energy of the diffusing photons. 
Then the typical energy of the photons decreases as they diffuse out. The typical temperature 
then falls below the photons temperature of the luminosity shell. If such processes exist at 
some radius r > r , how low can the temperature drop? Photons would only be created up 
until the radiation thermalizes. Given that the luminosity is constant, and that the diffusion 
time is rd/c, the photon energy density for all shells external to the luminosity shell up until 
optical depth of unity is 

e (1 < r < f) ^ (7) 
The thermal equilibrium temperature given such energy density is 

Tbb = (e/a)^/^ (8) 

where a is the radiation constant. The temperature will never drop below Tbb since then 
photons would be absorbed rather than emitted. As we show below, shells that can generate 
enough photons to affect the photons temperature on their diffusion out of the luminosity 
shell, will gain full thermal equilibrium, while shells that are out of thermal equilibrium do 
not affect the radiation temperature at all. 

To clarify, we provide an alternative description of the two options discussed above. 
As long as the radiation and the gas are in thermal equilibrium (i.e, e oc T^), the photon 
number flux increases with r. Thus, the gas at r > r must generate photons at a sufficient 
rate, and share the energy of the outgoing photons with the generated photons, in order to 
keep the radiation in thermal equilibrium temperature, Tbb- The observed temperature, is 
determined then by the outermost shell (at r > r ) that is in thermal equilibrium. If, on the 
other hand, none of the shells at r > Fis in thermal equilibrium then Tobs = T ^ Tbb- Note 
that we refer here to the typical photon energy as the radiation temperature also when the 
radiation is out of thermal equilibrium. 
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We note that the relation L = Airar'^T^, which is sometime used in the context of shock 
breakout, is wrong at any radius if kT is the typical photon energy at radius r {a is the 
Stefan-Boltzmann constant). The reason is that this equation, which relates luminosity, 
radius and temperature at the point where r ^ 1, assumes that the radiation is in thermal 
equilibrium at this point, while at the regime we discuss here it is never the case. At 
early times the optical depth is dominated by Thomson scattering, which conserves photon 
number. Therefore, as a typical photon should be absorbed at least once in shells that can 
keep thermal equilibrium, the shell at r = 1 is always out of thermal equilibrium and it plays 
no roll in determining neither the temperature nor the luminosity. 

Below we follow the photon-electron coupling in the different regimes and determine 
the time dependent observed temperature as function of initial conditions. We show that 
the observed temperature can assume radically different evolution between cases where the 
breakout shell is initially in equilibrium and cases where it is not. We stress that we use 
the term observed temperature to denote the typical energy of the observed photons also 
when the radiation is out of thermal equilibrium and the spectrum is not a blackbody. 
The observed temperature is often called color temperature when radiation is in thermal 
equilibrium (to separate it from the effective temperature, which we ignore throughout the 
paper since it is not an observable). We denote the shell at which the observed temperature 
is determined as the color shell. The color shell is the thermalization depth when thermal 
equilibrium is achieved, and as we show below the color shell is the luminosity shell when it 
is not. 



2.2.1. photon- electron coupling and thermal equilibrium 

The main processes that typically couple the gas and the radiation, at the physical 
conditions of interest, are free-free and bound- free emission and absorbtion and Compton or 
inverse Compton scattering. The importance of the free-free coupling, compared to bound- 
free coupling, depends on the metalicity, where in low metalicity environments the free-free 
process is more dominant while in high metalicity it is the bound-free. We assume here that 
free-free process dominate and discuss briefly the effects of bound-free couphng later. 

Consider an isolated shell with a radiation dominated energy density e, which generates 
its own photons by free-free emission. Now assume that initially, even though it is radiation 
dominated, it does not have enough photons to be in thermal equilibrium. The radiation 
temperature is T > Tbb and the hot photons must "cool"[^, by sharing their energy with 



^We use here the term "coohng" in the sense that the hot photons give their energy to cooler photons in 
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newly generated photons, for the system to approach thermal equilibrium. This is achieved 
by emission from electrons. The energy content of these electrons is negligible, but we 
assume that they can constantly gain energy from the photons an remain at the radiation 
temperature (this assumption is justified in appendix B). In this way, the energy of the 
existing photons is shared with that of the new photons that are constantly generated by free- 
free emission. Then, given a sufficiently long time, the number density of generated photons 
is hbb ~ CLT^^/SkTsB and thermal equilibrium is achieved. Therefore, the time required 
to achieve thermal equilibrium in the shell is roughly tibb I 'f^'phj fiTBs) ^ where Uphjf is the 
rate, per unit volume, at which free-free emission generates photons with energy hu ^ 3kT. 

In each shell that generates its own photons the time available to obtain thermal equi- 
librium is the time photons are confined to the shell, i.e., min{t, t^}. We, therefore, define a 
thermal coupling coefficient in the expanding gas: 

ubb 



r] = 



mm{t,td} fiphjfiTBs) 
7 ■ 10^ s / p V V ^Tbb ^ '^'^ 



(9) 



min{t,trf} V10-^°g/cmV \imeV 

where we approximate Uphjf ~ 3.5 ■ lO^^s^^ cm^'^ p^T^^/^. p is mass density (in g cm~^) 
and k is the Boltzmann constant. In the definition of rj we do not include Comptonization 
of low energy photons, which for the physical systems of interest becomes important only 
when the radiation is out of thermal equilibrium (see below), rj of the breakout shell at the 
initial time to is then approximated by 

/ \ 15/4 / N -1/8 

We took here into account the fact that the s hock compresses the gas by a fac tor of 7 and we 



used the relation Tbb^o ~ {Po'^o/'^Y^^ (see also lKatz. Budnik fc Waxmanll2010l ). The criterion 



for thermal equilibrium behind the sho ck at the breakout is ?7n < 1, which can be tran slated 



to a limit on the the breakout velocity (IWeaverlll976l : iKatz. Budnik &: Waxmanll2010l ): Vq < 



15, 000 km/s with a negligible dependence on po (as Po''^°). In case that bound- free emission 
dominates over free-free emission r^o can be smaller by up to an order of magnitude (see 
section 

For T] < 1, the radiation is in thermal equilibrium and the temperature is Tbb- Due to 
the relation between absorbtion and emission rj ^ 1 is also the requirement that a photon 
with hu = SkTBB is absorbed on average once by free-free process during the available time. 



the process. The total radiation energy is either constant, or reduced only by adiabatic expansion. 
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Thus, ii td < t then rj = 1 is equivalent to TabsT ~ 1 where Tabs is the free-free absorption 
optical depth for photons with hu = SkTsB- 

If 1] > 1 then free-free processes are not enough to couple the electrons to the radiation 
in the shell. Moreover, not enough photons are generated at hu ~ 3kT to obtain thermal 
equilibrium. However, if Compton scattering provides enough coupling so the electrons 
follow the radiation temperature, then the temperature T is determined by the total number 
of photons that are generated at hu ~ 3kT or that are generated at lower frequency, but 
can be Comptonized to an energy 3kT. Then free-free emission generates photons until 
the temperature satisfied e/3kT = mm{t,td} hphjf{T)^{T), where Uphj/ accounts only for 
generation of photons at hu ^ 3kT while ^(T) account for photons that are generated at 
lower frequencies and Comptonized to temperature T. Therefore, hphjj{T)^{T) is the total 
production rate of photons that can equally share the energy in the shell. Since free-free 
emission produces roughly a constant number of photons for every decade of energy below 
fcT, the Comptonization correction factor ^ is logarithmic. If no photons can be Comptonized 
then ^ = 1, while if Comptonization is important ^ may be much larger than unity. Writing 
this condition in terms of rj we find that if 77 > 1 then T^{TY = Tbbv"^- 



The availability of a low energy photon for computerization is discussed by IWeaver 



(119761 ). At the temperature and density range of interest (0.1 — 50 keV and ~ 10~ — 
10~^° gr/cm^) it is determined by the requirement that the photon energy can be doubled 
before it is re-absorbed by free-free proc ess. The rati o between T and the lowest energy 



photon that can be Comptonized is then (jWeaverlll976l ): 



_ kT _ ( p v^/y T 

^'""^ " - ^ ^10-9 g/cm-3 J \imev) ^^^^ 

Using the approximated rat e, per unit vol ume per unit of y = hu/kT, at which free- free 



emission generates photons (ISvenssodll984l ) we find that for Umax > 1 the number of gener- 
ated photons that can be Comptonized to kT is larger by a factor 1 (0.8 — ln[|/])— than 
the number of photons that free- free emission generates at ~ kT. Thus for Umax > 1: 

^(T) ^ max |l, iln[?/^aa;] (1-6 + ln[?/„a^])| . (12) 

If Umax < 1 then inverse compton does not contribute to the number of photons, implying 
^ = 1. All together the logarithmic correction to the number of photons, ^, ranges between 



•^Note that under the assumption that initially the energy density is radiation dominated, the system 
will always be driven towards this equilibrium since the "cooling" time of the system is shorter at higher 
temperatures (oc T~^/^). 
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^ = 1 at T < 100 eV to ^ ~ 100 at T = 50 keV. Therefore, since in the systems of interest 
Tbb ^ 100 eV Comptonization becomes important, and cannot be neglected, only once the 
radiation is out of thermal equilibrium. 

In summary, an isolated shell would arrive to a temperature of: 



1] > 1 



T = r^B^(^)^ ' (13) 

U 1] < 1 the spectrum is a blackbody. When rj > 1 the spectrum is an optically thin thermal 
free- free emission (with electrons at temperature T), which is modified by Comptonization. 
It is therefore a Wien spectrum (determined by the number of photons) at high frequencies 
and thermal spectrum (with T) at low frequencie^. 

Equation [T3] assumes that electrons and photons have the same temperature through 
the whole evolution. This is clearly the case after thermal equilibrium is obtained and while 
?7 < 1, since free- free emission and absorbtion keep the electrons and photons tightly coupled. 
In case that photons are out of thermal equilibrium they approach equilibrium by "cooling" 
via Compton scattering. This scattering is the heating source of the electrons, while the 
main electron cooling source is free-free emission. Therefore equation [13] assumes that the 
electron cooling is the slower process. In appendix B we show that this is generally true 
in the cases of interest, and that it may be marginally violated only in a small region of 
the relevant parameter phase space, in which case the modifications are likely to be small. 
Therefore we assume that electrons and photons have the same temperature in shells that 
generate their own photons and that equation [13] is always valid in these shells. 

Equation [13] neglects relativistic effects (pair prod uction, relatiyistic brernsstrah lung, 



etc.) and is, therefore, accurate only for T < 50 keV (jWeaverj [1976( ISvenssonl Il984[ l. We 
do not calculate here the exact temperature when eq. [I3] predicts higher temperature. In 
this case pair production will result in a lower temperature than eq. [13] predicts. In SNe, 
where the velocities are at most mildly-relativistic, the temperature should fall within the 
range of soft gamma-ray detect ors, ~ 50 — 200 keV, when pair production is important 



( IKatz. Budnik fc Waxmanll2010l ). 



The exact spectrum depends on the time it takes cold photons to double their energy by Comptonizition, 
(not calculated here). If photons have no time to be comptonized then the free- free emission spectrum is not 
modified. Otherwise a Wien spectrum, cx expl—hiy / kT], is observed in the region where the spectrum is 
dominated by inverse compton. Fi^ cx v'^ at lower frequencies, where the spectrum is dominated by free-free 
emission. At u < Vmm (defined in equation I lip the spectrum is dominated by self-absorbtion and F^ oa v^. 
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2.2.2. Succession of Shells & The observed photon energy 

In §2.2.11 we considered the evolution of the photon temperature in an isolated shell 
whose energy density e is known. We have assumed that the initial number of photons 
was small, so that the shell itself produced all its photons. We now discuss the interaction 
between shells as the photons diffuse out. A shell external to the luminosity shell receives 
photons from inner shells. It could, potentially, lower the typical photon energy by sharing 
the energy of the incoming photons with more photons that are produced in the shell. 

The effect of a shell on the radiation temperature depends on its thermal coupling 
coefficient: 

/ Tbb p \"Vmin{t,td}\"^ 

A shell cannot affect the observed temperature if it does not generate more photons than 
those arriving from inner layers. Since the energy flux, L, is independent of r at r > r, 
the number flux that each shell can generate is inversely proportional to the typical photons 
energy it could generate, given by equation [131 For 77 < 1 the number flux that a shell can 
generate is L/ikTBs) while for > 1 it is L/{kTBBif'i^'^)- Tbb is a decreasing function 
of r, and therefore shells with r/ < 1 increase the number flux while bringing the photons 
temperature to Tbb- On the other hand, as we show later, TbbV^ is an increasing function of 
r (the variation of the logarithmic factor ^ between shells can be neglected), and therefore, 
shells with rj > 1 cannot change the number flux of photons arriving from inner shells, and 
cannot modify the photon energies. Moreover, rj is an increasing function of r. Therefore, if 
the luminosity shell is out of thermal equilibrium, i.e., rj > 1, then all external shells have 
rj{m < m) > 1. It follows, that the only shell which affects the temperature when rj > 1 is 
the luminosity shell, which is generating its own photons. 

The observed temperature is therefore: 

T..^{i^:^ ^ r-' (15) 

[ TBBiV=^) V< 1 

Where ffc and are calculated for the shell that is currently the luminosity shell, at the 
point where the photon number in the shell was determined. Note that in case that rj < 1 
the color shell is exterior to the luminosity shell, where td < t. Therefore the criterion 
rj = 1 is equivalent to Tabs'^ ~ 1, which is a criterion for departure from thermal equilibrium 
used in stellar envelope c alculations, and in some of the works on shock breakouts (e.g.. 



Ensman fc Burrows! 119921 ). 



Tobs is the typical photon energy of the observed radiation, but the observed spectrum 
is not necessarily a blackbody. The spectrum of the radiation at the color shell is deter- 
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mined only by itself and is therefore that of an isolated shell with temperature Tobs (see the 
description below equation [13]) . As discussed above, outer shells do not affect the number 
of photons with energy hu > kTobs, however they may contribute to the number of photons 
with energy hp < kTobs, thereby modifying the low energy spectrum. We do not calculate 
here this modification. 

In the following sections we calculate the observed temperature by first finding Tbb in 
each shell. For shells external to the luminosity shell it is found using equation H] and the 
relation e ~ Lr/cr^. For inner shells it is determined by the adiabatic cooling during the 
expansion following the shock passage. Then we use equation [9] to calculate r] and equation 
[15] to find Tobs- We define the color radius to be the radius at which the observed temperature 
is determined (i.e., the radius of teh color shell), namely = max{r, r(r7 = 1)}. We define 
rrici as the color shell mass, i.e., m(r > r^). Note that we do not discuss the radius of the 
photosphere where r = 1, since it has no observable consequence or physical importance in 
the context of early SN light curves. 



As discussed above (section 12.11) . the optical depth of each shell is fixed during the 
planar evolution and the breakout shell is also the luminosity shell. Therefore, the observed 
temperature can only be determined by the breakout shell or by shells that are further out 
(i.e., m < mo). Unfortunately, the hydrodynamic evolution of these shells is uncertain, as it 
depends on the poorly understood details of the radiation dominated shock front as it breaks 
out of the envelope. As a working assumption, we take the density profile at m < ttiq during 
the whole planar phase to be similar to the pre-shock profile, i.e., p = pidQ/vot. We verify 
later that other reasonable profiles do not strongly affect the conclusions. Having p(m, t) we 
use equation [7] to find TBsifnyt) and then equation (noting that here td < t) to obtain: 



Thus, the thermal coupling increases slowly, but monotonically, with mass and time, ev- 
erywhere. This is true for any profile where the density drops with radius. The number 
of photons in the luminosity shell increases continuously if it is out of thermal equilibrium 
(i.e., fj decreases with time so % = fi{t)), and following equations [TSl and [T6] the observed 
temperature at to < ^ < is: 



2.2.3. The temperature during the planar phase (t^ <t < ts) 




(16) 




(17) 
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The value of Tobs,o when r^o > 1 is calculated by solving the implicit equation [13] for T using rjo 
and ^oiTobsfl)- Note that for > 1 equation [T7] is implicit due to the dependence of C, on Tobs- 
We solve for ^ by plugging p = po(^o/^) in equation [TT] ^ is continuously decreasing during 
the planar phase, but its decay is not described by a power-law and so is the decay of Tobs- 
Nevertheless the temperature instantaneous power-law decay index a = d\ia(Tobs) / d\ia(t) is 
bounded between that of an adiabatic cooling and that of a constant ^, i.e., 1/3 < a < 2/3. 
The drop in ^ is smaller when t] is larger, implying that a decreases with r]. 

The radius of the color shell is roughly constant during the planar phase, ~ -R*. If 
r/o < 1 then rrici depend strongly on the unknown density profile at m < niQ. If r/o > 1 the 
photospheric mass is constant rrici = wiq. We find that TbbV ^ increases with the radius, 
in agreement with our assertion that shells external to the breakout shell that are out of 
thermal equilibrium, do not affect T^bs since they generate less photons than those arriving 
from inner shells. Finally, we verify that taking other density profiles external to the breakout 
shell such as a steeper decreasing density (p = pidi/vot) and a constant (in radius) density 
(p = podo/vot) do not significantly change equation [T7] (the power-law indices when rf < 1 
are changed by less than 0.1). 



2.2.4- The temperature during the spherical phase t > tg 

The evolution of the spectrum during the early spherical phase depends strongly on 
whether the breakout shell is in thermal equilibrium at ts- When it is, the color shell at tg is 
at m < niQ and it quickly moves inward (in Lagrangian sense) to m > mg. Following similar 
steps to those we used to derive equation [TBI taking p(mo < m,t) we find 



-'-/^ / ^ \ 12(1.19n + l) / fJl \ 8(n + l) 



42n+49 22.32n+17 

r]{mo <m)=r]o., . , . 

We consider here only m < m where the time available for equilibrium is t^- The color shell 
is at m{ri = 1) and the observed temperature i^: 

8.03n+6 

2(1.76n+l) /+ \ 22.32n+17 



2(1.76n+l) /j- \ 



f \ (1.19n+l)(22.32n+17) 

(19) 



^ we neglect here the shghtly different time evolution between tg and t(riQ — 1), which introduces a 
negligible correction during this time. 
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The color radius and the mass above it are 



f \ (1.19n+l)(5.58n+4.25) 

r,i{r]<l)^R,{-] (20) 



and 

mci{r] < 1) = mo r]^ 



n + l 

2(i+l) f f \ 3(5.58u+4.25) 
5.58n+4.25 / ^ > 



7(6n + 7)(n + l) 
+ \ 6(1.19n+l)(5.58n+4.25) 

t) • 

The optical depth at the color shell decreases with time roughly as oc mlJrci oc 
^-0.4 _ Qj2ce Tci = 1 the opacity is not dominated anymore by Thomson scattering and our 
model assumptions break. Typically, however, the model assumptions break first by the 
temperature droping to 1 eV (and recombination becoming important) while > 1. Note 
that rci{ri < 1) increases faster than r with time. Thus, if the luminosity shell is at thermal 
equilibrium at a given time it will remain in equilibrium at any later time. Note also that 
equations [T9ll2T] are valid whenever 77 < 1, regardless of the value of rf at earlier times. A 
schematic light curve of this case is depicted in figure [H 

A very different evolution of To^s takes place if the luminosity shell is out of thermal 
equilibrium at tg- Here evolution at the beginning of the spherical phase is determined by the 
fact that the thermal coupling of each shell increases (i.e., 77 decreases) while its evolution is 
quasi-planar (i.e. t < R^/v). On the other hand the thermal coupling becomes weaker during 
the spherical phase. Therefore, the photon production of shells that are out of equilibrium 
(i.e., ?7 > 1) is negligible during the spherical phase and the number of photons in such shells 
is determined by rjc = rjmin, i-e. the minimal value of r] obtained when the shell just doubled 
its radius, and its thickness increased from the initial width di to -R* . There, using equation 
[HI with {TsB/TobsflY ~ {di/ R^Yf^mv^ / rriQVQ and p/po ~ fnto/mots, "we obtain (noting that 
during the planar evolution for the relevant shells t < td), 

I tn\ I m \ 6(n + l) 

Vmin = vit = Rjv)=r]oif\ i—j . (22) 

Thus, the observed temperature drops very quickly at the beginning of the spherical phase 
until ti = t(j]min = 1): 



ti — tf. 



+ \ 1/6' 

to 



3(1.19n + l) 
9.88n+5 



(23) 
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Between tg and ti the value of is also dropping significantly until it becomes of order unity 
at ti . Therefore Tobs does not evolve as a power-law during this time, and we do not provide 
here a formula for its evolution. Nevertheless the drop in T is very quick (oc t~^'^ when 
the evolution of is neglected) providing a strong observational signature, regardless of the 
exact decay rate. 

At ti the observed radiation is in thermal equilibrium (i.e.. Tots = Tbb), although rf > 1, 
since the radiation was in thermal equilibrium at t = R^:/v and equilibrium is kept through 
adiabatic cooling. The observed temperature during this phase, rfmin < 1 < vi't), is (note 
that = 1 at this time): 

-2/3 



+ \ ^'^ /f-,\ 3(1.19n + l) / / \ 6(1.19n + l) 

ToUh < t < h) = TsB,o vn -^] 7^ i^] (24) 



toy \ts J \ti 

12.48n+l 



The relation Tots = Tbb holds until ^2 = ^(^= !)• Since fj oc t6(i i9"+i'' we find 

1/6 



— t.Q 



6(1.19n+l) 
12.48n+l 

(25) 



During the whole time that rf > 1 the color shell is also the luminosity shell {vd = r) and 
therefore: 



0.81n + l 



and 



l \ 1.19n + l 

rci{ts<t<t2) = R,[j\ (26) 



2(n + l) 
^ \ 1.19n+l 

mci{ts < t < t2) = niQ (-j . (27) 



At t > ^2 the color shell is exterior to the luminosity shell (r < r^;) and is in thermal 
equihbrium. Thus Tobs, and rrici at t > t2 are given by equations (1911211 Finally, In all the 
different regimes of the spherical phase Tbbv"^ increases with the radius at r > r , implying 
that shells which are out of thermal equilibrium in this radii range do not affect Tobs- A 
schematic light curve in case that the luminosity shell is out of thermal equilibrium at tg is 
depicted in figure [TJ 



2.2.5. The effect of bound-free absorbtion and emission 



In our calculations above we assumed that free-free is the dominant emission and ab- 
sorbtion process. This is the case in low metalicity progenitor envelopes (^0.1 solar) but not 
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Fig. 1. — Main: A Schematic plot of the temperature evolution in the cases of ?7o > 1 (i^^d 
line) and ?7o > 1 (blue line). The power- law indices indicated in the plot depend weakly 
on n (see the text) and the values given here are very good approximations when n ranges 
between 1.5 and 3. Inset: A Schematic plot of the bolometric luminosity evolution. The 
luminosity, unlike the observed temperature, is independent of the thermal coupling and the 
value of 7]. 



in high metalicity ones (> solar). When bound-free dominates, the absorbtion coefficient, 
Kbf, is a non-trivial function of the wavelength, temperature and density. Nevertheless, a 
rough estimate of the opacity a.t u ^ kT/h show that the general deperidence on T and p 



is similar to that of free-free, i.e, Hbf oc hzjf oc pT~^-^ (e.g.. ISchwarzschildl 119581 ). Therefore, 
when bound-free dominates, thermalization is achieved faster (and kept for longer) by a 
factor Kbf/i^ff implying that in each shell 1] is reduced by this factor. Thus, all the equations 
derived above are valid after the definition of r] (eq. [9]) and the value of rjo (eq. [ID]) are 
divided by Hbf/i^ff (ignoring the effect of the logarithmic factor ^). 

The correction to the observed temperature in case that the radiation is in thermal 
equilibrium (r/ < 1) is small. The reason is that the dependence of Tobs on rjo is very weak 
(roughly as rj^^^), while in typical envelopes, which are dominated by Hydrogen or Helium, 
Kbf/nff < 10. This factor, however, is important in determining whether the radiation is in 
thermal equilibrium during the breakout and in calculating the temperature when it is not. 
Note that when the temperature is very high, > 10 keV, the metals are fully ionized and 
bound-free process can be ignored regardless of metalicity. 
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Fig. 2. — Main: Contours of 770 = 1 in the M^,-R^ plane for three different values of the 
explosion energy (3 ■ 10^" erg, 10^^ erg and 3 ■ 10^^ erg), and two values of n = 3 [solid 
line] and n = 1.5 [dashed line]). We used k = 0.34 cm^g~^ and assume that free-free process 
dominates the emission and absorbtion. Above and to the right of each contour, i.e. the 
larger and more massive progenitors, rjo < 1 and the radiation is at thermal equilibrium at 
the shock breakout time. Progenitors to the bottom and left of the contour (less massive and 
more compact) are out of thermal equilibrium during the breakout, and typically thermal 
equilibrium is gained only after the transition to the spherical phase. Typical M*-/2^, of 
RSG, BSG, WR and WD progenitors are marked as well. RSGs (n=1.5) are in thermal 
equilibrium during breakout of typical SNe and are expected to fall out of equilibrium only 
if the SN kinetic energy is > 5 ■ 10^^ erg. BSGs (n=3) are on the borderline and larger BSGs 
are expected to be in thermal equilibrium during breakout while less massive and compact 
BSGs don't. WRs (n=3) and WDs are out of thermal equilibrium in typical SNe breakout. 




3. The initial pulse 



In previous sections we calculated the luminosity and temperature of the expanding gas 
as a function of time. The initial timescale is to, and both the luminosity and the typical 
photons energy evolve as power- laws of time for t > Iq. However, several effects can smear the 
observed flux at early times. Differences in the light travel time to the observer, asphericity 
of the explosion (e.g., ICouch et al. II2OIOI ). or diffusion of the radiation through an optically 
thick surrounding such as thick wind blown by the progenitor before the explosion (e.g., iLi 



20071 ). Below we discuss how light travel time and stellar wind shape the luminosity and 



spectrum of the initial pulse. 
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3.1. Luminosity 

3.1.1. Light travel time of spherical explosion 

Light travel time in a splierical explosion significantly affects the light curve and spec- 
trum only at t < /2*/c and it would be important only if i?*/c > to- The light crossing time 
is always shorter than tg {R*/c < -R*/fo = tg), and therefore, its effects are important only 
during the planar phase. If light travel time shapes the initial pulse then for R^:/c > to, the 
observed flux rises over a timescale to and remains roughly constant over a duration R^^/c. 
The observed luminosity during the initial pulse is ~ Lotoc/R^ < Lq. After this plateau, 
at t > i?*/c, the luminosity starts decreasing as t~^^^ and is given again by equation as 
light travel time matters no more. An example of such light curve can be seen in figure [5] 
where the WR X-ray luminosity is bolometric until tg = 90 s. 

Note, that when light travel time shapes the initial pulse, then both i?* and to can be 
directly measured. The former is given by the duration of the initial pulse and the latter 
by its rise time. For more extended RSG progenitors, where to may be comparable to, or 
exceed, R^: (see appendix A) both the rise time and duration of the initial pulse are given 
by to and the progenitor's radius is not directly measurable. 

The discussion above is focused on the bolometric luminosity. For observations at 
frequencies lower than the initial typical energy given by Tobs,o, the rise time will be longer 
than to (assuming to < R*/c). Depending on thermal equilibrium it is either i?*/c or the 
time that the typical temperature falls into the the observed frequency window (see e.g., 
figures [3] and H]) . 



3.1.2. Winds 

WR progenitors are surrounded by the thick stellar wind ejected during the WR phase. 
The optical depth of a wind is gained mostly close to its source, between R^ and 2R^. Typical 
WR winds, are mildly optically thick, with optical depth that can be as hi gh as r,„ ~ 1 — 10 



once the wind is fully ionized by the precursor of the breakout emission (iLiI 120071 ) . If the 
wind is very thick, ^ c/vgh the radiation dominated shock would propagate in the wind 
as well and it should be treated in a similar way to our discussion in the previous sections. 
li 1 < < c/vsh then photons diffuse through the wind without generating a radiative 
dominated shock. The energy output of the shock breakout is not affected but the arrival 
time of the photons to the observer is smeared over T^R^jc. The pulse rise and decay times 
are both Tu,R^/ c. The information about to is lost and so is the ability to measure i?* directly. 
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3.2. Spectrum 

3.2.1. Light travel time 

If light travel time shapes the initial pulse, then at first, t to, the spectrum is domi- 
nated by the emission from the shock front which is propagating in a decreasing optical depth 
during the breakout. The typical observed photon frequency is therefore Tobs,o- As time 
evolves, to < t < R^^/c, the spectrum broadens to include lower frequencies as well. During 
this time high frequency breakout photons, Tobsfi, continues to arrive from areas with longer 
travel time while lower frequency photons from the expanding gas are arriving from areas 
with shorter travel time. Ignoring light travel time effects we derived Tobs = Tobs,o(t/to)~°', 
(1/3 < a < 2/3), while L = Lo(t/to)~^^^ (see equations (fT71) and (jlj)). Thus, the spectrum 
of the initial pulse broadens in time to form a power-law, 

oc (28) 

over a frequency range that grows with time. Its upper frequency corresponds to the initial 
temperature Tobs,o while the lower end of the power-law corresponds to the current (non 
delayed) temperature To{,s,o(V^o) The integrated spectrum of the initial pulse will show 
this power-law over a frequency range kTobs,o{cto/ R^:)°' < hu < kTobs,o- 



3.2.2. Wind 

In the case of a mildly opaque wind (1 < Tu, < cjv) photons spend a time tR^/c diffusing 
through the wind, thereby erasing all the temporal details on shorter time scales. As a result 
the observed spectrum is given by oc z/^^^ over the frequency range kTobsfl{cto / t^R*)'^ < 
hv < kTobsfl right from the beginning. 

Compton and inverse Compton scattering in the wind may also modify the photon's 
energy. However, the number of collisions per photon is r^, and the number of collisions 
needed to significantly change a photon energy is rrieC^/kT. Since the winds we are dealing 
with have moderate optical depth and our temperature calculations are applicable to cases 
where T < 50 keV, scattering within the wind cannot make a significant change to the energy 
of such photons. Therefore this effect is not very important in the temperature range that 
we can calculate. 
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4. Early SNe light curves from various progenitors 

Below we present early light curves (luminosity and spectrum) for different SN pro- 
genitors. We consider different progenitors of core-collapse SN: red supergiant (RSG), blue 
supergiant (BSG) and Wolf-Rayet (WR) stars. We also discuss the effect of deviation from 
thermal equilibrium on the sign al that follows the shock breakout of type la SN presented in 



Piro. Chang fc WeinbergI (120091 ). The light curves are derived according to the results pre- 



sented in previous sections, assuming that free-free is the dominant emission and absorbtion 
process, and the properties of the breakout shell given in appendix A. 



4.1. Red supergiant 

RSG is the progenitor of several members of the type II SN family. It has a convective 
envelope and its structure can be approximated using n = 1.5. We assume a hydrogen 
envelope with cosmic abundances so the scattering cross-section per unit of mass is k = 
0.34 cm^g~^ (the dependence on k is weak). We consider a typical radius of 500Rq (a light 
crossing time of about 20 min) and a typical mass of M^, = 15Mq. Following the initial pulse 
the luminosity evolves as: 



^^^^ ^ 3-10''eTg/sM^r'R^ooElf%'-'' t > U ^ ^ 



where = R^/xRq, = M^/xMq, = E/10^ erg and thr [td] is time in units of hours 
[days]. E here is the explosion energy, not to be confused with the internal energy in the 
expanding shells. The transition between the planar and spherical phases takes place around 



= 14 hr M^,,^^Rl-^E^r (30) 

The value of the thermal coupling parameter at the breakout is r/o = 0.06 M^^^''^'^ R^qq^^ E"^]^^ . 
Therefore, the observed temperature is determined at the outermost shell which is in thermal 
equilibriunj^ and it is given by eqs. [T7] and [T9j 

_ r 10 eV M,-,'-^'RliiEl,^Xr i<ts .... 



^In extreme cases (e.g., very energetic explosion with E^^i > 5), 770 may be larger than unity and the light 
curve would show a similar evolution to the one discussed in the context of BSG out of thermal equilibrium 
(see below). 
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Fig. 3. — The optical light curve (in absolute V magnitude) following the shock breakout 
from RSG [solid line), BSG in thermal equilibrium [thick dashed line) and out of thermal 
equilibrium during the breakout [thin dashed line) and WR [dash-doted line). In all cases 
the explosion energy is 10^^ erg. The progenitors radii are 500 (RSG), 70 R© (BSG 
thermal), 20 R0 (BSG non-thermal) and 5 R© (WR). The progenitors masses are 15 Mq 
in all cases except for the BSG thermal, where it is 25 Mq. The explosion is assumed to 
be spherical (light travel-time effects are included) and the pre-explosion stellar wind to be 
transparent. The source luminosity (before light travel time effects are included) at t < to 
is approximated as Loexp[l — y] (this is the luminosity of a radiation that leaks from the 
center of a static slab with a diffusion time to). The two thermal breakouts (RSG and BSG 
thermal) show a rising flux over a duration of R^/c. The rise does not stop at to since 
the optical band is below Tobs,o (the optical is not the bolometric luminosity at early time). 
At i?*/c the light curve start decaying until tg, when the flux start rising again during the 
spherical phase up to the time when Tots drops to the observed frequency (not seen here 
because we cut the hght curves at Tots = 10"^ K, see below). In the nonthermal breakouts the 
flux is strongly suppressed before tg since the temperature is very high (much higher than 
the optical). Then it rises very sharply up to ti, when it falls into thermal equilibrium. At 
t2 the flux joins the evolution of the thermal breakouts. Thus in the nonthermal breakouts 
the flux is rising continuously until the time that Tabs drops to the observed frequency. We 
present light curves up to the point that Tots = 10"^ K, since this is roughly the point where 
recombination, which is neglected in our model, becomes important. 




Fig. 4. — The FUV {i^p,jY = 2 • 10^^ Hz) light curve in absolute AB magnitude, following 
the shock breakout from RSG {solid line), BSG in thermal equilibrium [thick dashed line) 
and out of thermal equilibrium during the breakout {thin dashed line) and WR {dash-doted 
line). The parameters and assumption are the same as in figure [31 The light curves evolution 
is similar to optical ones, with the difference that the second peak in the flux is observed 
earlier, once Tobs get into the FUV range. 



The ratio of the diffusion time of the breakout layer and the star light-crossing time is 
cto/-R* = 0.25M°5^^i?|5{f E'^i^'™, implying that if the explosion is spherical the initial pulse 
has a rather well defined observed temperature of Tq^jisg ~ 25 eV M^^'^ R^qq^^ E^f and its 
rise time and duration are not too different. 

Typical optical and FUV light curves are depicted in figures [3] and |H The initial rise- 
time in both bands is i?*/c (assuming a spherical explosion). Note that unlike the bolometric 
luminosity, here there is no plateau between and -R*/c since the temperature is above the 
observed band at this time. Thus to cannot be easily recovered from optical/UV light curves. 
The temperature remains above both the optical and the FUV bands during the planar phase 
and both light curves decay slowly. During the spherical phase the luminosity drops more 
slowly while Tobs drops faster, as a result a break in the light curve is observed at t ~ and 
the optical fiux starts rising. It peaks once Tobs drops into the observed frequency range. The 
optical peak takes place after about two weeks. At that point the temperature is low enough 
so recombination, which we neglected, begins to be important. Therefore we terminate the 
light curve in the figures at earlier time when T^bs = 1 eV. 
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Fig. 5. — The total luminosity at frequencies > 0.2 keV {thick lines) and > 0.5 keV [this 
lines) following the shock breakout from RSG {solid lines), BSG in thermal equilibrium 
{dashed lines) and WR {dash-doted lines) stars. The parameters and assumption are the 
same as in figure [31 The depicted luminosity is the bolometric one, at t < -R*/c in the 
WR and BSG cases, showing the light curve of the initial pulse in the case of a spherical 
explosion. The luminosity rises until to when it becomes almost constant. The decay of the 
initial pulse starts at -R*/c. 



The X-ray light curve is depicted in figure O 7i?s'G(^o) is below the X-ray range and 
therefore only the initial pulse may be observed in soft X-rays. We assume here that the 
spectrum of the photons in the breakout layer is thermal, which implies that the initial X-ray 
pulse is very soft as the X-ray probes the exponential tail of the spectrum. Nevertheless, 
the energy in the soft X-ray flash from an RSG breakout is ~ 3 ■ 10^^ erg, which may be 
detectable to a substantial distance (possibly larger than that of a WR breakout X-ray flash). 

Note that our analysis is appropriate only for RSGs with a density profile that drops 
sharply with d near the edge, so the shock is accelerating before breakout. Our analysis 
also assumes that the width of the breakout shell is much sm aller than the stellar radius 
{dQ <^ R^). Stellar s tructure models suc h as those calculated bv lGezari et al. I (I20081), usiri g 
the KEPLER code (IWeaver et al.lll978l ) and 3D hydrodynamics code (IFreytag et al.ll2002l ). 
are compatible with these assumptions. Our model assumption s , how ever, are incompatible 
with density profiles such as the one used by ISchawinski et al. I (120081 ) . where the density of 
the envelope is very low (~ 10~^^gr/cm^ at 0.5i?*) and do ~ i?*. 
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4.2. Blue supergiant 

BSG is also a progenitor of member(s) in the type 11 SN Family. It has a radiative 
envelope and is well approximated using = 3. We assume a hydrogen envelope with k 
similar to RSG. The typical radius is SOi?© giving -R^,/c ~ 2 min. Following the initial pulse 
the luminosity evolves as: 



where tmin is time in minutes. The transition between the planar and spherical phases takes 
place around 

= 0.5 hr M^^RlfE^^^-^^ (33) 

The value of the thermal coupling parameter at the breakout is r^o ~ 2 M^^^'^"^ R^q'^ E'^f* . 
Therefore, BSG progenitors are on the thermal coupling borderline. Thus, we expect both 
types of BSG shock breakouts, where less energetic explosions of more extended and massive 
BSGs will be in thermal equilibrium, while more energetic explosions of compact less massive 
BSGs will be out of thermal equilibrium. For example an E = 10^^ erg explosion of a 
-R=K = 7O-R0 and = 25M0 progenitor has rjo ~ 0.5, and is therefore predicted to be in 
thermal equilibrium while if R^, = 20Rq and M^, = IOMq then 770 ~ 10 and the breakout is 
predicted to be out of thermal equilibrium. Finally, in high metalicity envelopes bound-free 
emission dominates over free-free emission and rjQ is reduced (section 12.2.51) . Therefore we 
provide here the temperature evolution for both thermal and non-thermal cases. 

If the radiation is in thermal equilibrium at breakout time, rjo < 1 (equations. [17] and 
[T9l) the observed temperature is: 



50 eV M,-,'-''R%^^El,^H-fJ'' t < t. 

f-0.11 pO.38 771O.IIJ--O. 
'25 -"-70 -^51 '^hr 



TBSG{Vo<l)-i lllyMl'-^'R^^^^li^H^'' uVt ■ ^^^^ 



In case that the radiation is out of thermal equilibrium at breakout time, 770 > 1, we 
use eqs. ([HD, ([21]) and ([ISD to find: 

Tbsg{vo > 1) 

r 150 eV Mro^i?2-o°-'^5¥Cn' t<ts 

2.1 



70 eV MfsO-^i^^oO-^Ei,! (±) ' < t < t, , (35) 

15 eV M?,^^R%^^E^r (1^)"° ' h < t < 



-27- 



where 



13 min M, 



--0.24 pO.94 77.0.29 
10 -'^20 -^51 



(36) 



and 



t2 = 20 mill M] 



— 0.77p0.62 rpO.99 
10 -"-20 ^51 



(37) 



as given by equations ( l23l) and (!25|) . The temperature at t < ti (when the observed radiation 
is out of thermal equihbrium) depends on the value of the logarithmic Comptonization factor, 
^, which does not evolve as a power-law. Thus, the temperature decay during the planar 
phase and early spherical phase {t < ti) is not an exact power-law. Here we provide a power- 
law approximation, using an average temporal indices a^^ and {TssGit < ts) oc t~"os^ 
Tsscits < t < ti) (X t"""!). The values we use here, a^^ = 0.39 and a^-^ = 2.1, are for 
the specific choice of the canonical explosion parameters, obtained by solving equation [T71 
numerically. As discussed above 1/3 < aos < 2/3, and is increasing when the breakout is 
farther from equilibrium. If for example we take an explosion of -E51 = 5 (keeping Miq = 1 
and R20 = 1) we find an average value of a^^ = 0.48 and = 2. Due to the dependence on 
^ the power-law dependence of TBscit < ^i) on the explosion parameters (Mio, -R20 and E^i) 
is also approximated and calculated numerically. The power law indices of Mio, -R20 and £"51 
that we provide are accurate to within ±0.5 in the range Miq = 0.3 — 3, R20 = 0.75 — 5 and 
E51 = 1 - 2. 

The ratio of the diffusion time of the breakout layer and the star light crossing time is 
cto/-R* = 0.09M°5^^i?55^^_E^j^°''^^. Thus, if r/o < 1 the observed temperature in the initial pulse 
smoothed between Tq bsg and 0.5Tq^bsg, while for r^o > 1 it is smoothed between Tq bsg and 
0.2To^BSG, where 



Similarly to equation |3ll The power law indices of Mio, -R20 and E^i in case that r/o > 1 are 
approximated and calculated numerically, with the same accuracy. 

Typical optical and FUV light curves are depicted in figures [3] and ID The light curves in 
case that thermal equilibrium is assumed are similar to the RSG case (with a different values 
of -R^k/c and tg) and they show an initial rise over a duration R^,/c, followed by a slow decay 
up to tg. In the non-thermal case the emission is strongly suppressed (as Tobs ^ Tbb) and the 
optical/FUV flux is rising sharply up to ti (slightly after tg) where Tohs = Tbb- As expected 
once Tobs = Tbb the two light curves, thermal and non-thermal, show a similar evolution. 
Note that the temperature at a given observer time during this phase depends mostly on 
the progenitor radius and therefore the peak in the light curve, which is observed once Tobs 
drop to the observed frequency, is observed at earlier time for more compact progenitor. 




m < 1 
Vo > 1 



(38) 
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The X-ray light curve when thermal equilibrium is assumed at all time, is depicted in 
figure [5l Tq bsg is ~ 0.1 keV and the X-ray luminosity is comparable to the bolometric 
luminosity. If the breakout emission is out of thermal equilibrium Tq^bsg is higher and the 
observed x-ray spectrum is harder. But in both cases most of the breakout luminosity, 
~ 3 • 10^^ erg, is expected to fall within the X-ray range. Making the X-ray flash from BSG 
breakouts the easiest to observe among the different progenitor types. 



4.3. Wolf-Rayet 

The core collapse of a WR star is most likely the onset of a type Ib/c SN. At WR 
stage the star radius is only several solar radii (i?*/c ~ 10s) and its mass is 10-80 Mq. 
The envelope is radiative (n=3) and has no hydrogen (k = 0.2 cm^g~^). In cases where 
the WR wind is very thick it may be dense enough to provide signiflcant opacity and play 
a role during the short planar phase (see section [3]). Note that this effect vanishes during 
the spherical phase, where r increases while the wind opacity decreases. Neglecting possible 
wind opacity effects the luminosity following the initial pulse evolves as: 



2 • 10^2 g^g/s 



0.33 d2.3 77'0.34j--4/3 f ^ r 
-"-5 -^51 ''min ^ f-s 



3.5 ■ 10^1 erg/s 



-0.73 



15 



(39) 



where 



90 s Mfs^i/^i-^SEslO-ss 



(40) 



The value of the thermal coupling parameter at the breakout is ?7o = 24 M^^^'^"^ R^^'-^ E'^f^ , 
implying that initially the radiation is not in thermal equilibrium with the gas and the 
observed temperature is given by eqs. f|T7|) . flMl) and f|T9|) : 



WR 



-2.2 



( 1 keV Mr/-^/?5-°-2E5Vt-°-' 
140 eV Mfs^-^i^o-^Eii^ 

40 eV Mri?r^^5l°"Cn' 
5 eV Mr5°-iii?0-38E0iiHr°-6i 



^51 '^hr 



t < ts 

ts <t <ti 

ti <t <t2 
t2 < t 



(41) 



where 



and 



200 s Mr5°-''/?g-''i?5°f 



400 s Mr5°-''i?g-^'i?5°r 



(42) 



(43) 



Similarly to the nonthermal BSG case the temperature at t < ti does not evolve as a power- 
law, so the values a^^ = 0.4 and = 2.2, are approximations obtained for the speciflc 
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choice of the canonical explosion parameters (see the discussion below equation [37]). For the 
same reason the power-law dependence on the explosion parameters at t < ti is approximate, 
where the power-law indices of M15 and are accurate to within ±0.5 and of E^i to within 
±0.8 in the range M15 = 0.2 - 2, i?5 = 0.2 - 2 and E51 = 1 - 2. 

The observed temperature at the first few minutes is much higher than the prediction of 
a model that assumes thermal equilibrium at all time, peaking at the breakout temperature: 

To,wR ~ 2 keV Mr,'-' R^' '^ Elf (44) 

Similarly to Tobsit < ts) the dependence of Tobsfl on the explosion parameters is approximated 
and calculated numerically with a similar accuracy to the indices provided for Tobs{t < tg). 
Finally, Note that our model, which neglects relativistic effect, is not applicable at T > 
50 keV. Thus whenever the model predicts higher temperature it overestimates the true 
temperature which does not exceed ~ 200 keV. 

In a typical WR star cto/R^ = 0.014M°527i?o.9i^-o.73_ Therefore the initial pulse is 
always composed of a range of temperatures (regardless of the wind opacity). The integrated 
spectrum shows a non-thermal power-law, F,^ oc z/~^ where < /3 < 0.5, over more than one 
order of magnitude at hu < kTo^wR- 

Typical optical and FUV light curves are depicted in figures [3] and IH The flux in these 
bands is very faint at t < ts (as long as the emitted radiation is out of thermal equilibrium) 
and it rises continuously until Tabs drops into the observed band. Figure [5] depict the observed 
flux in X-rays and soft gamma-rays. This energy range contains almost all of the breakout 
luminosity. The radiation from less energetic breakouts from larger WR progenitors will 
all be in the range of X-ray detectors (0.2-10 keV), while more energetic breakouts from 
compact WR progenitors will be in the range of soft gamma-ray detectors(> 10 keV). Thus, 
we expect that at least some WR breakouts can be detected as soft gamma-ray bursts by 
satellites such as Swift up to a distance of ~ 10 Mpc. 



4.4. White dwarf 



The thermo-nuclear explosion of a white dwarf near the Chandrasekhar mass is most 
likely the o rigin of a type la SN. The sho ck breakout in this case was carefully explored 
recently by iPiro. Chang fc Weinberg! (120091 ). This paper assumes a thermal equilibrium at 
all time, while pointing ou t that this assumption may be v iolated at early times. Here we 
repeat the main results of iPiro. Chang fc Weinberg] (120091 ). highlighting the effects of the 
deviation from thermal equilibrium at early times. 
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The shock velocity at breakout may be either mildly relat i vistic {vq/ c[1 — Vq/ c^]~'^^'^ ^ 1) 
or relativistic ([1 — Vq/c^]'^/'^ ^ 1) ( Tan. Matzner &: McKeell200ll ). Here we consider only 
cases where the breakout is mildly relativistic and therefore a Newtonian approximation of 
the hydrodynamics is reasonable for the breakout properties, in such case tg ~ -R*/c and 
the energy in the initial pulse is Eq ~ 10^*^ — 10^^ erg over a duration R^/c ^ tg ~ 10 — 20 
ms. The initial diffusion time is much shorter (to ~ /isec), implying that if asphericity 
does not play a major role, the bolometric luminosity rises practically instantaneously. The 
temporal evolution of the lu minosity at t > t,, is independen t of the assumption of thermal 
equilibrium and is given in iPiro. Chang fc WeinbergI (120091 ). The normalizations that we 
provide here, based on the simple assumptions described in appendix A (which ignores the 
phase during which the nuclear burning provides energy t o the shock), are comparable to 
the more accurate calculation of iPiro. Chang fc WeinbergI (120091 ): 



LwD ~ 

5 . W^'-fM^;i 
3 ■ 10392S 



0.73 R 
5-108 cm 



0.35 



/\^-0.73 R 77-0.91 ( tdeg \ 
^^Ha 5-108 cm-^51 yiOOs ) 



-0.18 



, tg <^ t <^ t(^e 
^-0.17 . f ^ ^ 



(45) 



where t^eg is the time in which m makes the transition from mass that was not degenerate 
be fore the explosion (n = 3) to ma tter that was degenerate (n=1.5). In the model discussed 



bv lPiro. Chang fc WeinbergI ( 120091 ) tdeg ~ 700 s. 



The breakout temperature, assuming thermal equilibrium is Tbb,o ~ 10 keV, but rjo ~ 
10^ implying first that the system is far from thermal equilibrium and second that relativistic 
effects, such as pair production, play a major role. Therefore we cannot determine the exact 
observ ed spectrum at early time, but p air production limits the temperature at this regime 
to be ( IKatz. Budnik fc Waxmanlboiol ): 



Tf 



0,WD 



200 keV 



(46) 



This is the temperature in the shocked gas frame, and it is also the observed temperature if 
the shock is not relativistic. In case that the shock is relativistic the observed temperature 
is higher. Therefore a breakout from a type la SN will produce a short flash of gamma-rays 
which is easily detectable within our Galaxy by various gamma-ray satellites such as IPN 
satellites and Fermi, and may possibly be detectable out to the Magellanic Clouds by Swift. 
At t > ts the temperature drops quickly until it gets to Tbb at t2 ~ 1 s. The short flash of 
7-rays may be classified as a short GRB without the following detection of a SN. At later 
times the observed temperature is: 



TwD = 6 eV Ml 



-0.11 



R 



0.38 



5-108 



cm 



77.0.11^.-0.61 

51 ^min 



leg- 



(47) 
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This equation can be extrapolated to t > tdeg using Two oc t~^'^^. Note that only at t ^ 1 s 
the velocity of m is significantly lower than c (e.g., v{ls) ~ 0.5c), limiting the applicability 
of our temperature evolution calculation to t ^ 1. At this point the emitted radiation is 
already in thermal equilibrium. Yet, our predicted temperature at late times, whe n ther- 
mal equilibrium holds, does not agree with this of iPiro. Chang fc Weinbergl (120091 ). They 
assumed that the temperature is related to the luminosity and the radius of the luminosity 
shell, Tobs = {L/Anar^f/^. However, as we have shown, the photons that are emitted from 
the luminosity shell are reprocessed within all the external shells that are in equilibrium. 
Thus, we use Tobs = T(r] = 1) (see the discus s ion in section [2^2]) . and predict lower temper- 
atures compared to lPiro. Chang fc Weinbergl ( 120091 ). Eq. HT] implies that the light curve we 
provide is valid only during the first ~ 20 min. At later times the temperature is low enough 
so recombination can no longer be neglected. 



5. Comparison to previous works 

Various aspects of the early supernova light curve were already addressed by a large 
number of authors, both analytically and numerically. The early SNe light curve that we 
consider here is composed of three phases: breakout, planar phase and spherical phase. 
Almost all analytical calculations focused only on a single phase. Therefore we first compare 
our methods and results to previous analytical calculations phase by phase. Later we compare 
our results to those of numerical calculations that include (implicitly) all the phases and the 
transitions between them. 



5.1. analytic works 



The method we use to find the breako ut shell, i.e., finding Tq = c/vn, is similar to the 
one u sed in previous analytical studies (e.g.. llmshennik fc Nadezhirull988uMatzner &: McKee 
19991 ) and so are the values of Eq and Lq that we find. Also the temperature of the breakout 
shell is similar, but only if the breakout shell is in thermal equilibrium. The observed 
breakout temperature depends on whether the post shock radiation in the breakout shell is 
in thermal equilibrium or not. If it is, then the color shell is exterior to the breakout shell but 
interior to the layer where r = 1. Therefore, the observed breakout temperature, Tobs,o, is 
lower than Tbb,o but larger than the effective ter aperature, defined as T^ff = \L/4:Tiar J_-iY/^. 
This fact was recognized by several authors (e.g.. lFalklll978l : lEnsman &: Burrowslll992l ). They 
estimate numerically that during breakout the ratio of the observed temperature, (denoted 
sometime as the color temperature Tc), to the effective temperature is Tobsfl/Tefjfl ~ 2. 



- 32 - 



Similarly, in our analytic theory we obtain that this ratio is 1.8 for our canonical RSG and 
2.1 for our canonical BSC 

We determine if the breakout shell is in thermal equilibrium when the shock breaks out, 
by equating the time available for photon production, min{t, t^}, and the time needed to 
obtain thermal equili brium (we use the same method also at any other shell and time). This 
method was used by IWeaverl (119761 ) to analytically find departure from thermal equilibrium 
just behind radiative dominated shocks, where the available time is roughly the shock crossing 
time. During breakout t ~ in the breakout shell and both are co mparable to the shock 
crossing time. Therefore, our estimates should coincide with those of IWeaverl (119761 ) when 
the properties of the shock once it breaks out are taken. Our requirement for departure from 
thermal equilibrium at breakout, r^o > is translated at the densities of interest (the density 
dependence is extremely weak) to energy per nucl eon > 1 MeV or a shock velocity > 15, OOP 
km/s. This result is similar to the one obtained by lWeaverl (119761 ) (his equation 5.15). [Weaver 
(Il976l ) also calculated numeric ally the temperature in case tha t the gas behind the shock 
velocity is out of equilibrium. iKatz. Budnik fc Waxmaru (120101 ) realized the importance of 
this result to SN shock breakouts and presented an analytic calculation (their equation 
18), which is similar to our equation [T3I (the two equations are almost identical when the 
numerical factor in our equation [10] is taken as 0.4 instead of 0.2 and our ^ is identified as 
their Aeff qeff). Our analytic calculation is in excellent agreement with the numerical results of 
Weaver! (119761 ) over the whole range of relevant densities and velocities before pair production 
becomes important. 



The only analytical calculation of the planar phase was carried out by lPiro. Chang fc Weinberg 



(120091 ). in the context of SN type la breakout from a wh ite dwarf. The luminosity we ob- 
tain is similar to that of jPiro. Chang fc Weinberg J2009h . Howe ver, the observed temper- 
ature that we find is different. jPiro. Chang fc Weinberg] (120091 ) assumes thermal equilib- 
rium and Toft,, = Tbr- At ear ly time thermal equilibrium is not achieved and therefore 



Piro. Chang &: Weinberg (l2009l ) underestimate the observed temperature. At late time the 
radiation thermalized further out than the luminosity shell and therefore they overestimate 
the observed temperature. A detailed comparison with their work is presented in §4.4[ 

The spherical phase was explored analytically, using different methods by a number of 



authors ( 


Chevalier 


1992; 


C] 


levalier & Fransson 


2008; 


Waxman, Meszaros & Campana 


2007; 


Rabinak & Waxman 


20ld 


Piro. Chang & Weinberg 


2009 


). We calculate the luminosity by 



identifying the energy source of the observed radiation as the point where p hotons can diffuse 
out over a dynamical time . This method use the same physical picture as I Chevalier! (119921 ) 
and I Chevalier fc Fransson] (120081 ) . which carried out a more detailed calculation by finding 
self similar solutions for the diffusion wave that propagates into the ejecta. Their results are 
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similar to ours. IWaxman. Meszaros fc Campanall2007l and lRabinak fc Waxmaru (120101 ) use a 
different method to calculate the luminosity. They use the equation L = Anar'^^^T^^^ (which 
assumes thermal equilibrium at the r = 1 shell), where T^=i is calculated by pretending that 
the radiation in the r = 1 shell has cooled adiabatically since breakout. However, the 
radiation that was generated by the shock at that shell is long gone by the time that r = 1 
(it escaped once its opacity satisfied t = c/v) and the radiation that arrives to the r = 1 shell 
from inner shells is not in thermal equilibrium. This method overestimates the luminosity 
by a factor of r ^-o^-r for ri = 1.5 — 3. For the specific problem and parameters that we 
consider here r = 30 — 70 during the spherical phase of a core-collapse SN, translating to an 
overestimate of the luminosity by up to a factor of 2. 

All previous works assumes that the observed radiation is in thermal equilibrium during 
the spherical phase. We find that thermal equilibriu m is always achieved long time after the 
transi tion between the planar and spherical phases . IChevalierl ( 1l992l ): IChevalier fc Fransson 
( 120081 ) and IWaxman. Meszaros &: Campanal (120071 ) assumes that the radiation is in thermal 
equilibrium all the way out to the location where r = 1. This leads to an underestimate of 
the observed temperature by a factor of r(?7 = 1)*^'^^, which for our canonical SN parameters 
is an underestimate by up to a factor of 2 during the spherical phase. 



Rabinak &: WaxmanI (120101 ) explore SN light curves at T < 3 eV, including the effect 



of recombination and taking into account Thomson, free-free and bound-free opacity. Our 
calculation, which ignores recombination, can be compared to t heir treatment only at T > 1 
eV, where the deviation due to recombination is less prominent. iRabinak &: WaxmanI (120 lOf ) 
look into the difference between the observed temperature and effective temperature. They 
find Tobs/Teff{l < T < 3 eV) ~ 1.2, which is slightly lower than the values that we obtain 
(1.3-1.5) when bound-free opacity is ignored. 



5.2. numerical works 



Many authors used numerical simulations to study early SN light curves (e.g..lShigeyama. Nomoto &: H 



1988 



Wooslev 



1988 



Ensman fc Burrows 1992; Blinnikov et al. Ill998: Schawinski et al. 11200 



Tominaga et al. II2009I ). Two of these works provide bolometric luminosity and observed tem- 
perature at temporal resolution that is high enough for comparison wit h our model, starting 
at the breakout through the planar and s pherical phases. These are lEnsman fc Burrows 
(119921 ) that simulate BSG explosions and iTominaga et al. I J2009h ;hat simulate an RSG 
explosion. 



10^ 



10^ 
time [s] 



Fig. 6. — Bolometric luminosity {upper panel) and observed (color) temperature {lower 
pa nel) of two BSG exp l osions as a function of time, as obtained by the numerical simulation 
of lEnsman fc Burrows! (119921 ) {dashed black line) and by our analytic calculation {solid blue 
line). The progenitor parameters are M15 = 16/15 and R50 = 0.9. The two explosions 
energies are E51 = 1 {thin lines) and E^i = 2.3 {thick lines). The luminosity and temperature 
of the ■E'.Pii = 1 explosion are divided by a factor of 4 (both ours and those obtained by 
Ensman fc Burrowslll992l ). The temperature of the E^i = 2.3 explosion in case that thermal 



equilibrium is artificially enforced is also plotted {thin dashed- dot line). See text for details. 



Ensman fc Burrowsl (119921 ) simulate several BSG explosions. They provide detailed lu- 



minosity and temperature curves, with no corrections for light travel time effects, of a BSG 
progenitor with = 16 Mq and R^, = 45 R0 that explodes with energies of 10^^ erg 
and 2.3 ■ 10^^ erg. The temporal resolution is about 10 s for the first few hundred seconds, 
adequate for comparison with BSG breakout and planar phas e evolution. At later times 
they provide curves with a temporal resolution of about 0.5 hr. lEnsman fc Burrowsl (119921 ) 
allow for deviation of the gas temperature from the radiation temperatures and identify the 
color shell (i.e., thermalization depth) by the requirement y/3rabsT 
impose erad 



2/3. They do however 

where Crad and Trad are the radiation energy density and temperature 
respectively, thereby not allowing for deviation of the radiation from thermal equilibrium. 
Figure [6] sh ows the bolometric l u minos ity and observed temperature as functions of time 
obtained by lEnsman &: Burrowsl (1l992l ). These are extracted from their figures 2 & 11 at 
early times and 3 & 12 at late times, t = is set such that the peak of the luminosity 
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is at to. Figure [6] also shows our results of the luminosity and observed temperature. Our 
luminosity is calculated using equation 139) for both explosion energies. The excellent agree- 
ment between Ensman fc Burrows Jl992f )" numerical luminosity calculation and our analytic 
results is achieved without any fitting of the normalization or any other parameter. The 
results range over three orders of magnitude in time and luminosity and do not differ more 
than a factor of 2. The temporal decay slopes are similar and so are the values during the 
planar and spherical phases. 

For E^i = 1, the breakout thermal equilibrium is marginal (r/o ~ 2) and we plot only the 
ter nperature curve expected i f thermal equilibrium is kept at all time (as enforce artificially 
by lEnsman &: Burrows! Il992l ) using equation [311 For E^i = 2.3, the breakout radiation 
is out of thermal equilibrium {rjo ^ 13) and therefore we calc ulate the temperatur e using 
equation [351 In order to compare our model to the results of [Ensman &: Burrowsl (Il992[ ) 
we also plot the temperature predicted in that case if there were thermal equilibrium at all 
time, using equation [Ml T he agreement between our model (when thermal equilibrium is 
enforced) and the results of [Ensman fc Burrowsl (Il992l ) is very good. It is better than 20% 
at any time, except for a period around a few hundred seconds where the temperatures of 



Ensman fc Burrowsl (119921 ) drop over a short period by a factor of two. We do not see any 



physical origin for this behavior, which may be a numerical artifact. In the more energetic 
explosion, E^i = 2. 3, our model shows that t he assumption of thermal equilibrium fails at 
early time and that [Ensman fc Burrowsl (119921 ) underestimate the breakout and planar phase 
observed temperature by an order of magnitude. 



Tominaga et al. I (120091 ) simulate a 1.2 ■ 10^^ erg explosion of an RSG progenitor with 
= 18 Mq and i?* = 800 R0. They provide bolometric luminosity and observed tem- 



perature curves, where a travel time effects are included, assuming a spherical explosion. 
The temporal resolution is about 0.5 hr for the first 14 hour and about 1 day at later 
times. We extract these curves from their figure 3 and set t = so the luminosity peaks 
at t = _R,,/c = 1866 s. Figure [7[ sho ws the bolometric l umin osity and observed tempera- 
ture as functions of time obtained by [Tominaga et al. \ (120091 ) and those obtained by our 
calculations. Our luminosity is calculated by using equation [291 and smoothing it for photon 
arrival time effect from a spherical explosi on. The result is then di vided by a constant factor 
of 1.4 to best fit the luminosity found by [Tominaga et al. I (l2009l ). The similarity between 
the result of the detailed numerical simulation and our calculation is very good. Again, the 
luminosity ranges over three orders of magnitude and it is similar to within a factor of two 
at all times (factor of 3 if our luminosity is not divided by a constant factor of 1.4). The 
observed temperature is calculated using equation [HH At early times, when light travel time 
effects cause the observer to see a range of temperatures at any given moment, the maximal 
of these temperatures is taken. Our resulting temperature is not multiplied by any constant 
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Fig. 7. — Bolometric luminosity {upper panel) and observed (color) ter nperature (lower 



pa/net) of an RSG explosion as function of time as obtained numerically by lTominaga et al. 



( I2OO9I ) [dashed black line) and by our analytic calculation [solid blue line). The explosion 
parameters are M15 = 1.2 and i?5oo = 1-6 and E^i = 1.2. The luminosity calculated by our 
analytic formula (equation [29]) is deviled by a constant factor of 1.4. see text for details 



factor. The agreement between the numerical and analytic temperatures is again very good. 
It is similar at any time to within 30%. 



6. Summary 

We derive analytic SNe light curves at early times, as long as recombination and radioac- 
tive decay do not play an important role. These light curves are valid while the observed 
temperature is above about 1 eV and before injection by radioactive decay becomes impor- 
tant. These conditions hold during the first day after the explosion of a typical SN. The 
main advantage of our analysis over previous ones is the account for the radiation-gas cou- 
pling, which leads to determination of the observed temperature when the radiation at the 
color shell is out of thermal equilibrium. It also corrects previous estimates of the observed 
temperature, and the color shell location, when the radiation at the color shell is in ther- 
mal equilibrium. We define a thermal coupling coefficient, rj, and find that the temperature 
evolution can follow two very different tracks, depending on rjo, i.e., the value of rj in the 
breakout shell at the breakout time. When the breakout shell is out of thermal equilibrium 
(?7o > 1) the observed temperature starts high above the value obtained when thermal equi- 
librium is assumed and it drops faster than the case that the breakout shell is in thermal 
equilibrium. Thermal equilibrium is typically gained (when rjo > 1) only at early stages of 
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the spherical phase. 

We discuss the luminosity and spetral evolution during the initial pulse and derive early 
SN light curves for various SN progenitors as a function of the explosion energy and the 
progenitor mass and radius. These are useful for interpretation of SNe light curves during 
the first day, which can teach us about properties of the progenitor star and potentially to 
lead to its identification. Additionally, it can be used to evaluate the effect of the early 
emission (e.g., ionization of the circum burst medium), in case that its detection is missed, 
on the environment at the SN vicinity. Finally it is useful for planning targeted searches of 
shock breakouts of various SNe types. The theory we discuss here can be also applied in 
some cases to shock breakout from non-SN stellar explosions. For example, the explosion of 
sol ar-like star by tida l force s in the vicinity of a super-massive black hole discussed recently 
bv lCuillochon et al. I J20091 ). 



The main conclusions based on our analysis are: 



It was shown that shoc k breakout radiation from WDs, WRs and some BSGs is out 
of thermal equilibrium ( iKatz. Budnik fc WaxmanI 120101 ) . We show that it typically 
remains out of thermal equilibrium throughout the planar phase and until the early 
spherical phase. In SN from these compact progenitors the observed temperature at 
this time is significantly higher than the one obtained when thermal equilibrium is 
assumed. The observed temperature falls as where 1/3 < a < 2/3, during the 
planar phase, and once the evolution becomes spherical it plunges down (roughly as 
until we observe radiation that is in thermal equilibrium at the source. 

Breakouts from RSGs and some BSGs are in thermal equilibrium. The flux at frequen- 
cies below Tohs (e.g., optical/UV) starts with a bright initial pulse and then it decays 
during the planar phase reaching a minimum at tg. The flux is rising during the planar 
phase reaching a second maximum when Tobs falls into the observed frequency. 

In cases where the radiation is in thermal equili brium at the source, the l ocation of 
the thermalization depth, is not trivial (e.g., lEnsman fc Burrowsl Il992l ). The as- 
sumptions used in some pr evious analytic c alculations, such as Vri = rjr = 1) (e.g.. 



Chevalier &: Franssonll2008l ) or r^i = r (e.g., 
rect. Instead, in this case < < r{T = 



Piro. Chang fc Weinberg! |2009| ) are incor 



1) and it satisfles Vd = r{r] = 1), where 



The initial pulse, which is smeared by light travel-time and asphericity, contains emis- 
sion from gas at different temperatures. The resulting combined spectrum is not a 
black-body. Instead, below the peak of uFj^ the spectrum is a power-law oc u'^ 
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where < /3 < 1/2. If asphericity is negligible then this effect is negligible in RSGs, 
minor in BSGs and significant in WRs, where a power- law oc z/~^, ranges over one 
or two orders of magnitude in v. 

The bolometric luminosity is dictated by the balance between the radiation diffusion 



time scale and the dynamical time scale (e.g.. IChevalierlll992[ ). Namely, the source of 
the observed energy at any time is the shell where these two timescales are compa- 
rable. Since the photons dominate the heat capacity, the bolometric luminosity does 
not depend on the photon-gas coupling, and it evolves independently of the thermal 
coupling. 

The bolometric light curve and spectral evolution of the initial pulse depends on the 
explosion asphericity and wind transparency. A "bare" spherical explosion has a unique 
and therefore identifiable light curve. It shows multiple time scales as well as spectral 
evolution during the initial pulse. It rises quickly, over t^, and lasts for R^/c. A 
mildly opaque wind {1 < < c/vq) results in a pulse of characteristic time Ty^R^/c 
during which the spectrum does not evolve significantly. A "bare" aspherical explosion 
produces a time evolving spectrum, whose hght curve and spectral evolution depends 
on the asphericity details. Much can be learned from observation of the initial pulse! 

The initial pulse of a RSG may release ~ 3 ■ 10^^ erg over a duration of ~ 1000 sec in 
soft X-rays ~ 0.2 keV. A BSC shock breakout releases a comparable amount of energy 
in harder X-rays over ~ 100 sec. The breakout from a WR releases ~ 10^^ erg within 
~ 10 sec in the form of hard x-rays or soft 7-rays. Thus, not only WR stars, but 
also breakouts from BSGs and potentially RSGs are predicted to produce strong X-ray 
flares that can be detected by X-ray telescopes to large distances. WR breakouts may 
also be hard enough to be detectable by current soft gamma-ray detectors to a distance 
of ~ 10 Mpc, and mimic a single pulse long gamma-ray burst. 

The early emission from a type la SNe is in 7-rays. It is detectable within our galaxy 
by current detectors and may mimic a very short (10 — 20 ms) gamma-ray burst. 
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APPENDIX A- parameters of the breakout shell 

Here we give a short list of values of various physical parameters of the breakout shell 
at the time of breakout. We consider three different progenitors of core collapse SNe - RSG 
(n=1.5, K = 0.34 cm^g-i), BSG (n=3, k = 0.34 cm^g-i) and WR (n=3, k = 0.2 cm^g-i). 

We approximate the density profile near the star surface as: 

. (A-l) 

where = M^/B?^. There is a correction factor of order unity to eq. lA -!( which de- 
pends on the stellar struct ure, but the results are practically insensitive to this factor 



( ICalzavara fc Matznerl 120041 ). The properties at different locations near the stellar surface 



are 



m = l^r^) (A-2) 



n + 1 VP* 

r = ^f^) (A-3) 



n + l 



\ 1/2 / N -0.19 
^51 \ 



Note that we ignore the difference between the ejected mass (used ,e.g., by lMatzner fc McKee 



19991 ) and the total stellar mass. The following properties of the breakout shell at the time 



of breakout are found by requiring r = c/fQ. 



10-3 Mo M0543i?2^26^-0.56 (j^gg) 

mo^ { 3 ■ 10-6 Mq M0541/?2^34^~0.58 (BgQ) (A -5) 

3 ■ 10-8 Mq Mi^l ^2.34^-0.58 (^j^) 

6 ■ 10-10 g/cm3 M{'567i?-lg-64^-0.3 ^j^gQ) 

2 ■ 10-9 g/cm3 M{'556i?5-l-25^-0-44 (ggg) (A -6) 

6 ■ 10-8 g/cm3 M{'556i?-l-25^-0.44 (^j^) 

7000 km/s M^^-^^'R^^^'^Elf (RSG) 

20, 000 km/s M{^-^^Rl^-^^Elf (BSG) (A -7) 

40, 000 km/s M'^-^^ Rf-'^'^ E^f (WR) 

9 ■ 10^^ erg M^^^-^^Rl^^Elf (RSG) 

E^^ { 3 • 10^6 erg M^^-^^Rl^^Elf"" (BSG) (A -8) 

9 ■ 10^^ erg M^^-'^^Rl-^'^Elf (WR) 
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Note that Eq is the internal energy in the breakout shell at to, while E^i is the total explosion 



energy in units of 10^^ erg 



50 MiV3^0o26^-0.56 (j^gG) 

ro^<; 15 M?5^ii?of ^5-^8 (BgG) 
8 MOg^ii^o.ss^-o.ss (^j^) 



(A -9) 



^0 



5.8 • 10-3 M,t''RtolE^r 



6.1 • 10- 



-0.15 



1.9 ■ 10-3 M-°-i^i?0-58E-°-i^ 



^51 



(RSG) 
(BSG) 
(WR) 



300 s M^,^^Rl,l^E^r 

10 s MOg^^i^l^Sl^-O.TS 

0.2 s MOg^^i^i.gi^-o.rs 

0.25 M°,^ii?^o^o^E5-i°-'' 
0.09 M^s^^i^agi^-o-rs 
0.014 MOg^^i^o.gi^-o.rs 

14 hr M^,fRl,l^E^r 

0.58 



0.5 hr Mri?^o''^5i 



90 s MiV'i^^-'^^^sl"-'' 

0.06 Mr/-^^i?,-oy^E,v^ 
2 Mr5^-6^i?,-oi-iE,Y^ 



24 Mi 



-1.63 



-1.1 



7-12.24 
^51 



obs.O 



25 eV Mr5°-^i?5-oo''^5V 
80 eV M25°-2«i?7-o°-''^5V' 
700 eV Mfo^-^i^-'-iEsV 



2 keV Mfs^-^ 



-1.5 



-'^51 



(RSG) 
(BSG) 
(WR) 

(RSG) 
(BSG) 
(WR) 

(RSG) 
(BSG) 
(WR) 

(RSG) 
(BSG) 
(WR) 

(RSG) 
(BSG thermal) 
(BSG nonthermal) 
(WR) 



(A -10) 



(A -11) 



(A -12) 



(A -13) 



(A -14) 



(A -15) 



Note that Tobs,o is the observed temperature at the time of breakout and not necessarily the 
temperature of the breakout shell at that time. The dependence of Tobs,o in the nonthermal 
cases (BSG nonthermal and WR) on M*, R^: and E^i is approximated and calculated 
numerically. The reason is that the value of the logarithmic Comptonization factor ^ has a 
non-power-law dependence on the explosion parameters. The power law indices we provide 
are accurate to within ±0.5 when i?* and are larger or smaller by up to a factor of about 
3, than our canonical values and E^i = 1 — 2. In order to obtain a more accurate evaluation 
of nonthermal Tobsfl when M*, i?* and E^i differ from our canonical values, 7po (the factor 
of 7 is due to the shock compression) should be plugged into equation [11] iymax) and then 
equation [T3] should be solved numerically (using r]^). 
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APPENDIX B - Electron-photon temperature coupling 

Our temperature calculations (in particular equation [T3l) assume that in the relevant 
shells the electron temperature follows the photon temperature at all time. Here we examine 
the validity of this assumption. 

Consider first a shell that is in thermal equilibrium at the time of breakout (i.e., ?7j < 1). 
As discussed in section 12.2.11 when 77 < 1 a typical photon is absorbed at least once within 
the available time {min{t,td})- Therefore, since photons dominate the energy density, the 
electron temperature is coupled to the radiation while 77 < 1. Now, if in a shell rji < 1 then 
by the time that rj > 1 that shell does not affect the observed temperature anymore since 
at this point it is external to the luminosity shell. Therefore in such shells our temperature 
coupling coupling holds while relevant. 

Next consider a shell that is initially out of thermal equilibrium (i.e., rji > 1). Here, 
absorbtion do not play a significant role and electrons follow the photons temperature if 
the rate at which they lose energy (electron cooling rate), mostly via free-free emission, 
is slower than the rate at which they gain energy via Compton scattering with typical 
photons (electron heating rate). The coupling in shells that are out of thermal equilibrium 
is important mostly during the planar phase, since during the spherical phase the color shell 
is quickly receding inward to shells that are in thermal equilibrium. Since during the planar 
phase only the breakout shell is observed (it is also the luminosity shell and the color shell 
when f] > 1) we consider only the coupling in this shell. Below we calculate the heating 
and cooling rates in the breakout shell that is out of thermal equilibrium during the planar 
phase. 



The electron heating term is: 

cuph SkT 
I rip rrtpC 



eueat = i^^.^kT (B-1) 



where here we denote number densities as n with different subscripts, not to be confused 
with the stellar envelope structure power-law index, e is the energy of a single electron, 
/ = l/rieCT is the photons mean free pass, Uph is the photon density, rie is the electron 
density and we take 3fcT photons as responsible for most of the heating. Just behind the 

2 

shock Uph^i ~ where n„ is the upstream density (compressed by a factor of 7 in the 

shock, i.e., Ue^i = Tn^), implying: 

eheat ~ caTUe—^—kT (B-2) 

where E is the total energy in the radiation and Ei is its initial value. During the planar 
phase E/Ei oc f-^/^. 
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The main cooling process is free-free emission: 

ecooi ~ aarcue \/m^c^kf (B-3) 

where a is the fine structure constant and we neglect the logarithmic factor ^. We define 
the heating to cooling ratio: 

X - ^ - 4 (^^) ' (^) (L\ (B-4) 

If X > 1 then electrons track the photons temperature, while if x < 1 then Compton 
scattering is not enough to keep the electrons at the photons temperature. Thus, If x < 1 
and in addition r] > 1 m the breakout shell during the planar phase, then not Compton 
scattering nor free-free absorbtion can heat the electrons fast enough and our assumption 
of a single temperature fails. Using equations [HI [16] and kTsBfl ~ ('^oPo/'^)^^'^ we find that 
while 7] > 1 then in the breakout shell: 

Thus, our assumption of a single electron-photon temperature is valid during the breakout. 
Moreover, x depends strongly on rj and only weakly on t (e.g., in a typical BSG, (tg/toy^^ ^ 
8). Therefore, the request x > 1 while rj > 1 may be at most marginally violated in the 
breakout shell during the planar phase, and only in a narrow range of the parameter phase 
space. Note that even if the electron temperature in the breakout shell decouples from that 
of the photons during the planar phase the modifications to the observed light curve are 
minor. The bolometric luminosity is of course not affected while the temperature in the 
breakout shell is not driven anymore towards thermal equilibrium and it starts dropping 
only via adiabatic cooling, e.g.. Tots oc t~^^^ (instead of a slightly faster decay). Later, the 
temperature evolution during the early spherical phase is also slightly modified. Nevertheless, 
the general behavior, that the temperature drops very sharply at tg < t until it approaches 
Tbb and merges with the thermal equilibrium evolution at t ^ t2, is unchanged. 

Therefore, we conclude that our assumption that the electrons and photons temper- 
atures are similar generally holds in the shells that determine the temperature evolution, 
and that even if it is violated over a small region of the relevant parameter phase space the 
modifications are small. 



APPENDIX C - Glossary of main symbols and notations 

t: time since breakout 
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■ r: radius 

• v: velocity 

• m{r): mass at radius larger than r 

• p: mass density 

• d: shell width 

■ r: optical depth 

■ E: internal energy (not to be confused with E^i) 

■ td- diffusion time 

• e: energy density 

■ rj: Thermal coupling coefficient. Defined in equation [9] 

■ ^ : logarithmic Comptonization term. Defined in equation [12] 

• L: observed luminosity 

• Tobs- observed temperature (often denoted as color temperature). Defined as the typical 
photon energy. 

• breakout shell: The shell where the shock breaks out (the shock width is comparable 
to the shell width). 

• luminosity shell: The shell that generates the observed luminosity. In this shell the 
diffusion time equals to the dynamical time, t. 

■ color shell: The shell where the observed temperature is determined. This shell coin- 
cides with the thermalization depth when the observed radiation is in thermal equilib- 
rium and with the luminosity shell when it is not. 

• For any quantity x, we use the following subscripts and superscripts 
xf. Initial value (after shock crossing) of a shell. 

xq] value at the breakout shell at the time of breakout. 
x: value at the luminosity shell. 
Xd'- value at the color shell. 

■ to- Duration of shock breakout. Also, dynamical time, diffusion time, and shock 
crossing time of the breakout shell at the time of breakout. 
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■ ti: The first time where the observed temperature is in thermal equilibrium if it was 
not in thermal equilibrium from the beginning. Thermal equilibrium was achieved by 
photons produced at earlier time. Defined in equation [23] 

■ t2'- The first time where the observed temperature is in thermal equilibrium, which 
is achieved by photons produced at t2, if it was not in thermal equilibrium from the 
beginning. Defined in equation [25] 

• tg-. The transition time from planar to spherical geometry (= R^^/vq). 

■ Tobsfl'- observed temperature at to- Equal to the temperature of the breakout shell if it is 
out of thermal equilibrium. Otherwise it is lower than the breakout shell temperature. 

■ Tbb'- Thermal equilibrium temperature appropriate for a given energy density by Boltz- 
mann's law. Defined in equation [8] 

• i?*: stelar radius 

■ M*: stelar mass 

• E^i'. explosion energy in units of 10^^ erg 

• n: power-law index describing the pre-explosion stellar density profile near the edge. 
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